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Abstract 

Beyond point mass effects various contributions add to the radiative evolution of compact bi- 
naries. We present all the terms up to the second post-Newtonian order contributing to the rate 
of increase of gravitational wave frequency and the number of gravitational wave cycles left until 
I^"i the final coalescence for binary systems with spin, mass quadrupole and magnetic dipole moments, 

f^ ' moving on circular orbit. We evaluate these contributions for some famous or typical compact bina- 

f^ , ries and show that the terms representing the self interaction of individual spins, given for the first 

^SJ . time here, are commensurable with the proper spin-spin contributions for the recently discovered 

» , ' double pulsar J0737-3039. 
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1 Introduction 



Neutron star and black hole binary systems are among the significant sources of gravitational radiation 
which detection is expected by the gravitational wave observatories. The frequency of the radiation 
Aq ' emitted by these binaries is expected to be in the sensitivity range of the Earth-based interferometric 

fSJ I detectors p. Recently observations are under way to directly detect such signals and analysis methods 

^D ' were specified for inspiral signals from binaries with 3-20 solar masses and for setting upper limits on 

'^ I inspiral event rates for binary neutron stars using interferometer data P] ■ The detection of gravitational 

radiation from these compact binary systems is also expected by the Laser Interferometer Space Antenna 
f— ^ ! Q] ■ Alternative theories of gravity can be tested and parameters of spinning compact binaries can 

"^ ' be estimated from these measurements. 

^^1 The final coalescence of compact binaries is preceded by an inspiral phase for which the post- 

I ' Newtonian (PN) approximation provides a reliable description. This description is considered valid 

{^1 until the system reaches the innermost stable circular orbit. In |S] the authors have discussed the prob- 

lem of the failure of the PN expansion during the last stages of inspiral, called the intermediate binary 
black hole problem |7], both for spinning and non-spinning black hole binaries on quasicircular orbits. 
j^ ' Reliable results can be achieved by stopping the integration at the minimum of the energy as function 

H I of orbital frequency [H] . Tidal torques become important in latter stages of the inspiral E^ • 

The equations of motion were given to 3.5 PN order accuracy in the post-Newtonian regime with the 
inclusion of spin-orbit (SO) effects and their first PN correction in JT]. Spin-spin (SS) ^21i quadrupole- 
monopole (QM) ^31 and magnetic dipole - magnetic dipole (DD) contributions ^^ to the accelerations 
were also discussed. The rate of the radiative change of otherwise conserved quantities, like the secular 
energy and angular momentum losses {dE/dt) and (dJ/dt), characterizes the backreaction of gravita- 
tional waves escaping compact binaries on the orbit, which were computed to leading order in [151 116) 
and its PN and 2PN corrections in ^7] and |18| . 

When the interaction of the spins Si with the orbit is taken into account, the total angular momentum 
J = L + Si + S2 is still conserved, however the orbital angular momentum L is not, due to spin 
precession. Its magnitude L remains constant as a consequence of the specific functional form of the 
precession equation ^H]- The radiative change of (dE/dt) and (dL/dt) characterize the backreaction on 
the radial part of the motion which were computed in |19l [SOj . 

The SS contributions to (dE/dt) and (dL/dt) were derived in |^|^, where a description in terms of 
the magnitude of an angular averaged orbital angular momentum L was introduced. This method was 
used in the discussion of the QM [23| and DD contributions [23 . The magnitude L of the orbital angular 
momentum is not conserved due to the spin precessions caused by spin-spin, quadrupole-monopole and 
magnetic dipole - magnetic dipole interactions. 

The SS contributions in (dE/dt) and (dL/dt) given in [21] contained not only interaction terms 
between the two spins, but self-interaction spin (SS-self) terms. These are originated from the terms 



proportional to JgQ (bln) J^q (ajv) in dE/dt and e^-'^JgQ {aL]\j) Jg^ (ajv) in dJ^/dt, where Jg^ 
denotes the n*^ derivative of the spin-orbit contribution of the velocity quadrupole moment evaluated 
with the Newtonian acceleration a^r . 

Here we compute the SS-self contributions to the rate of increase of the gravitational wave frequency 
/ and to the accumulated number of gravitational wave cycles Af. For completeness we enlist all other 
contributions to df/dt and TV to 2PN order, namely the PN, SO, SS, QM, DD, 2PN and tail terms. 

Due to the emission of gravitational waves the orbit of the binary tends to circularize 1151 . Accord- 
ingly we consider circular orbits for which the gravitational wave frequency is twice the orbital frequency 
P5] . In the next section we evaluate the rate of increase of / which is given by the rate of change of the 
orbital angular frequency to — nf under radiation reaction: 



dt J \duj J \ dt 

where the expression dE/du is found by differentiating E = E{uj). The various contributions to the 
secular energy loss 

dE\_/dE\ /dE\ /^^\ /^^\ 

dt / \ dt / ^ \ dt / pjv \ dt / goj^tail \ dt / 2PN + (SS-self)+SiS2+QM+DD 

were given in C7]-|23 and P^ . 

Eq. (Q is immediately integrated since it is an ordinary differential equation in u. By a second 
integration we obtain the accumulated number of gravitational wave cycles Af. We evaluate the different 
contributions to Af for some famous or typical compact binary systems and compare the magnitude of 
the SS-self term with other contributions. 

2 Frequency evolution 

Introducing spherical coordinates the components of the acceleration a in the system (n. A, Ln) are 

n ■ a. = r — ruj , A • a = rtj + 2fu! , Ln ■ a = —ruj A • — — , (3) 

where r = rn is the separation vector, Ln = /ir x r is the Newtonian orbital angular momentum, 
A = Ln X n and a hat denotes a unit vector. The orbital angular velocity uj is introduced by the relation 
V = rn + rujX. For circular orbits r = r = and v^ = muj holds. Here m = mi + m2 is the total and 
/Lt = mim2/m is the reduced mass of the system. Consequently (mw) is of PN order. The radial 
projection of the acceleration yields the orbital angular velocity as r • a = —r^uP' (PJ. From the explicit 
form of the acceleration, with various contributions given in |12l 1131 I14j we find lo — lo (r) and then 
r = r (w), 

3 3 ^-^ \ m^ / mf 
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In Eq.|01) we have introduced the following notations. The magnitude of the spins and magnetic dipole 
moments are denoted by Si and di. In a coordinate system with the axes (c,L x c,L), where c is the 
unit vector in the J x L direction the polar angles of the spins are Ki and ipi. In [23 other coordinate 
systems were introduced with the axes (bi. Si x bi. Si) with bi are unit vectors in the Si x L directions, 
respectively. In this system the polar angles of the the magnetic dipole moments di are a^ and (3i. 
Moreover 7 — cos~-^(Si -82), A ~ cos~^(di • d2) and 77 — ^/m. The parameters Pi — Qi/mim^ 
characterize the quadrupolar contributions with the quadrupole- moment scalar Qi [I'd] . In the last term 
Ao is defined as 

^0 = 2cosA + 3(picr2 -p20-i)sin('(/'2 - V'l) - 3(piP2 + cricr2)cos(?A2 - V'l) , (5) 

Pi — sin ai cos /3i , (Ji = cos ai sin Hi + sin a^ sin (3i cos Ki . 



E = E (uj) is obtained by the combination of Eq. Q with the expression of the velocity on circular 
orbits V — ruj in the energy integral E = E [r,v): 



^H = -iM™c.)^/4l-J(3 + |)(mc.f/^ + mc.|]Q;;| + 2r; 



^(co) = -^/^(mw)^/^<;i-^(3 + ^)(mw)"/^ + mw>J( r^ + 2r;) ^cosk, (6) 



+ (mco)4/3 



-(-27+ 19?7 - — -) H T- (cos7 - 3cosKiCOSK2) + - > p.t(3cos Kj - 1) H — - 

8 3 rym* 2 ^^ 77m'* 



The radiative evolution of the orbital angular frequency for circular orbits is deduced from Eqs. Q, Q 
and lO as 



dt / ^ 5 
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where 

O" = f S1S2 + f^SS-self + CFQM + (^DD ■ (8) 

Here /?, csiSai '^SS-seif, ctqm and (Tdd are the spin-orbit, spin-spin, self-interaction spin, quadrupole- 
monopole and magnetic dipole-dipole parameters: 
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i=l 



^ = ^^4(ll3^ + 7577)cos^. , (9) 



SiS; 



'^SiS2 = -^^^(-247cos7-h721cosKicosK2) , (10) 

'^SS-seif = ^E(^) (6 + sin^-O , (11) 

2 — 1 ^ ^ 

5 ^ 
o'QAf = --^Pi (3cos^Ki - l) , (12) 

1=1 
5 
o-_D_D = jdid2Ao ■ (13) 

The N, PN, SO, SS, 2PN and tail contributions to Eq. ((TJ were verified to agree with the results of |77| . 
the QM and DD terms with J3I and |T3], respectively. 

Eq. (0 is an ordinary differential equation in u;. To linear order in the perturbations all angular 
variables can be considered constants. Since all the angles appear in perturbative corrections they 
are given with sufficient accuracy to PN order, in which order they do not change jl9|. Hence the 
time-evolution of uj for circular orbits can be obtained by an integration over time: 

.it) = ^-^llUJ^+'-l,]r-y^ + ^(^-.]r-y^ (14) 

^ ^ 8to [ V2688 32 7 10 ^4 / ^ ' 

I 1855099 56975 371 ^ 3a\ _i/2l 

''' 1^14450688 ^ 258048'^ "*" 2048 '^ ^ 64 y J' 

where the dimensionless time variable r = r]{tc — t)/5m is related to the time (ic — t) left until the final 
coalescence. The accumulated orbital phase 0c — f/* is given by a further integration: 



/ ujMdr 

V J 
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Table 1: The accumulated number of gravitational wave cycles. The frequencies fin and fjm are given 
for the frequency domain of the detectors or up to the innermost stable circular orbit. The first two 
columns refer to the LIGO/VIRGO type detectors and the last two to the LISA bandwidth. 



The accumulated number of gravitational wave cycles emerges as 

1 f ,/« /3715 55 



Af = 




8064 96 7 4 14 



/ 9275495 284875 1855 . 15cr\ ,,^ 
\ 14450688 258048 ' 2048 ' 64 J 



(16) 



The 2PN and tail contributions agree with those given in 28 and the other terms with Eq. (4.16) of 

We enlist all contributions to TV in Tablen]in terms of /3 and a, evaluated from Eq. 116|l for compact 
binaries. These are the double pulsar J0737-3039 29' I30L I5T] . one neutron star - stellar mass black hole 
binary and two examples of galactic black hole binaries [H] . 



3 Conclusions 

Various contributions add to the rate of increase of the gravitational wave frequency and the accumulated 
number of gravitational wave cycles. Here we have presented the complete set of contributions to the 
evolution of gravitational wave frequency and to the accumulated number of gravitational wave cycles 
up to the 2PN order, namely the PN, SO, SS, QM, DD, tail and 2PN terms, with the inclusion of the 
previously unknown self-interaction spin terms. These results add to the closed system of first order 
differential equations governing the secular evolution of radiating compact binaries already derived in 
|19[ I21L 1221 I23L 1^ and represent an important step towards a complete characterization of the orbital 
evolution. 

To comment on the importance of the self-interaction spin contributions we evaluate the spin pa- 
rameters CTSiSj and ass-seif for the double pulsar J0737-3039. The neutron stars in this double pulsar 
have average radii 15 km, masses of 1.337M0 and 1.25Af0, pulse periods of 22.7 ms and 2773.5 ms. By 
the Jenet-Ransom model |22] the angle ki (in their notation A) has two possible values: ki — 167° ± 10° 
(Solution 1), and ki = 90° ± 10° (Solution 2). The other angle K2 can be determined by solving nu- 
merically the spin precession equations which was done for black hole binaries in |33) . According to 
|34[ I35j it is likely that wind-torques from the energetically dominant component have driven the spin 
axis of the other component to align with the direction of L, causing K2 — 0. The estimates for the spin 
parameters are given in Tabled 

We conclude that the proper spin-spin contribution vanishes whenever one of the spins is aligned 
with the orbital angular momentum and the other spin is perpendicular to it. Thus crss-seif is the only 
spin-spin contribution in this case. Moreover, though four orders of magnitude smaller than the spin- 
orbit effects, already considered in "SE^ , self-interaction spin contributions are found to be comparable 
with the proper spin-spin contributions. This is due to the fact that one of the spins is two orders of 
magnitude larger than the other. 
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Spin parameter (order) Solution 1 Solution 2 

crss-seif (10 '^) 

Table 2: Spin parameters for the two solutions of the Jenet-Ransom model representing the binary 
pulsar J0737-3039 

The self-interaction spin contributions are more important when one of the spins is negligible com- 
pared to the other. The motion of a test particle orbiting around a massive spinning body is described 
by the Lense-Thirring approximation. Its first correction in the gravitational radiation is represented 
by the SS-self contribution of the higher spin. 
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